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Abstract 

The dependence of the dimensionless Brans-Dicke (BD) coupling parameter upon time and 
the scalar field has been determined, for an isotropic and homogeneous space of zero 
curvature, by solving BD field equations and the wave equation for the scalar field. For this 
purpose, very simple expressions of empirical scale factors, that generate constant 
deceleration parameters, have been used here in two theoretical models. The characteristics of 
time dependence of the BD parameter, obtained from these two models, are in qualitative 
agreement with each other. The mathematical expressions representing the rate of change of 
the BD parameter with time, based on these two models, are found to have identical forms. 
Combining the expressions of the BD parameter obtained from these two models, a method 
of determination of the present value of the equation-of-state (EoS) parameter of the cosmic 
fluid has been discussed. Its value is found to be consistent with the ranges of values obtained 
from other studies based on recent observations. 
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1. Introduction 

The information regarding several interesting phenomena in the field of cosmology, 
particularly the accelerating expansion of the universe, have emerged from astrophysical 
observations [1, 2]. The recent observations reveal that the accelerated expansion of the 
universe is caused by an exotic form of energy with negative pressure, referred to as dark 
energy (DE). This DE is one of the most important subjects of scientific investigation in 
physics and astronomy. To explain the late-time acceleration of the universe, lots of models 
have been proposed. The cosmological constant (A) is the simplest representative of DE [3]. 
But the models involving A are known to have cosmological constant (CC) problem (the fine- 
tuning problem) and the problem of cosmic coincidence [4], To explain the accelerated 
expansion by solving such problems, many alternative models have been proposed either by 
modifying the right hand side of the Einstein field equations by incorporating specific forms 
of energy-momentum tensor or changing the left hand side of Einstein field equations. The 
first category includes the DE models of quintessence [5], k-esssence [6], the family of 
chaplygin gas [7], holographic [8], new agegraphic [9] dark energy etc. The second category 
models are based on modified gravity, such as f(R) gravity, scalar-tensor theories [10] etc. 
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The Brans-Dicke (BD) theory of gravity is a path breaking study of scalar-tensor theory by 
Brans and Dicke, incorporating Mach’s principle into gravity [11]. In BD theory, the space- 
time dynamics is represented by the metric tensor and the dynamics of gravity is represented 
by a scalar field. The gravitational constant in BD theory is replaced by the reciprocal of a 
time dependent scalar field (<p). There is a dimensionless parameter (co) that couples this 
scalar field (< p ) with gravity. BD theory has successfully explained the experimental findings 
from the solar system [12], Brans-Dicke theory has been found to be very effective regarding 
the recent study of cosmic acceleration [13]. From general theory of relativity (GR), this can 
be obtained by letting oo -> oo and (p = constant [14]. The early and the late time behaviour 
of the universe have been very effectively solved by this theory. According to a study by 
Banerjee and Pavon, BD theory can potentially solve the quintessence problem [131. The 
generalized BD theory is an extension of the original BD theory where the BD parameter oo is 
regarded as a function of the scalar field <p [15]. The generalized BD theory gives rise to a 
decelerating radiation model where the big-bang nucleosynthesis scenario is not adversely 
affected [13]. A self-interacting potential have also been introduced in its modified form. 
Bertolami and Martins have used this theory to explain accelerated cosmic expansion for a 
spatially flat model [16], According to these theories oo should have a low negative value, of 
the order of unity, to account for the accelerated expansion. This is in contradiction with the 
solar system experimental bound of oo > 500. Bertolami and Martins have successfully 
explained the accelerated expansion with a potential (p 2 and large |co| [16]. 

In the present study, the time dependence of the BD parameter (m) has been determined from 
the field equations of the Brans-Dicke theory of gravity and the wave equation for the scalar 
field (<p). Here we have formulated two theoretical models of the BD parameter, based on 
two empirical expressions of the scale factor. One of them has a power-law type of 
dependence upon time and the other one has an exponential dependence upon time. Each of 
these very simple scale factors produces a time independent deceleration parameter. Although 
the deceleration parameter is known to be a function of time, these scale factors have been 
chosen for the simplicity of mathematical calculations, without much loss of generality. The 
scalar field, in each model, has been assumed to have a power law type of dependence upon 
the corresponding scale factor. The values of the constants involved in all these empirical 
relations have been evaluated from the field equations. The properties of time variation of the 
BD parameter, obtained from these two models, are found to be in qualitative agreement with 
each other. They are consistent with the findings of other resent studies in this regard [17, 
18]. Apart from time dependence, the dependence of the BD parameter upon the scale factor 
and the scalar field has also been determined in the present study. Although the entire 
formulation is based upon a time independent equation-of-state parameter, an attempt has 
been made, at the end of this analysis, to determine the present value of the EoS parameter. 
Combining the results obtained from the two models, a mathematical exercise has been 
carried out for an approximate estimation of its present value. The findings of this theoretical 
method are in agreement with the range of possible values of this parameter, obtained from 
recent astrophysical observations [19, 20]. This result seems to provide justification for the 
choice of scale factors that generate time independent deceleration parameters. 

2. Solution of Field Equations 

In the generalized Brans-Dicke theory, the action is expressed as [11], 

S = f d^x^g ^ gd^cp d v cp + (01) 
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Here, g is the determinant of the tensor metric g /iV , R is the Ricci scalar, L m is the Lagrangian 
of matter, (p is the Brans-Dicke scalar field and m is a dimensionless Brans-Dicke parameter. 
In generalized BD theory, this coupling parameter oo is considered to be a function of (p. 

The line element for a non-flat universe, in the framework of the Friedmann-Robertson- 
Walker (FRW) cosmology, can be expressed in polar coordinates as, 


ds 2 — —dt 2 + 


a 2 (t) 


dr + r 2 dd 2 + r 2 sin 2 d d% 2 

.1 — KT^ 


( 02 ) 


Here a(t) is the scale factor of the universe, t represents the cosmic time, r is the radial 
component, k is the curvature parameter and 9 and / are the two polar coordinates. 

The gravitational field equations of Brans-Dicke theory, for a universe filled with a perfect 
fluid and described by FRW space-time with scale factor a(t) and spatial curvature k, are 
given by, 
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2 a a 2 +k cjjcp) <p 2 

a a 2 2 (f> 2 


+ 2 + — 

a cp <p 


YP_ 

<P 


(04) 


The wave equation for the scalar field (<p), in generalized Brans-Dicke theory of gravity, 
where m is a time dependent parameter, is expressed as, 

ip ^cup _ p (1-3 y) <P <0 

<p acp <p 2CJ+3 <p 2o) + 3 V 


Combining equations (3) and (4), and taking k — 0 (for flat space), one obtains, 
2- + /^ — +S-- + -= p( - 1 ~ y - > 

a a 2 a cp cp (p 


(06) 


The energy conservation equation for the cosmic fluid is written as. 


p + 3^(l + y)p = 0 (07) 

In the above expressions, y (= P/p) is the equation of state (EoS) parameter for the cosmic 
fluid, with specific values for different cosmological periods. The values of y are, 1 (massless 
scalar field dominated era), 1/3 (radiation dominated era), 0 (matter dominated era), — 1 
(vacuum energy dominated era). 

Considering this EoS parameter y to be time independent, the solution of equation (7) is 
obtained as, 


p = p 0 a 3 ( 1+y ) 


(08) 


Combining equation (3) with (8) one gets. 


-,« 2 | Ppg 3(1+y) 

. . _ a 2 a<P _ <P 

1 ip* 

2 (p 2 

Combining equations (5) and (8) one gets, 


( 09 ) 
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1 ppa 3 ( 1+ ^ ( 1 —3y)—cLxp _ ' 

2 ij)+3d(p/a 


( 10 ) 


To determine the time dependence of a) from the above equations we have used two sets of 
empirical expressions for the scale factor (a) and the scalar field (cp) in the following two 
models. In these two models, the scale factors are such that one obtains constant deceleration 
parameters from them. The scalar field in each model has been chosen to have a power-law 
dependence on the corresponding scale factor. 

The values of different cosmological parameters used in this article are: 

72— 

H 0 = = 2.33 x 10“ 18 sec^, q 0 = -0.55, p 0 = 2.83 x l(T 27 tfg m~ 3 

(Po — — — 1.498 x 10 lo Kg s 2 m~ 3 , t 0 = 1.4 x 10 10 Years = 4.42 x 10 17 s 

G o 

MODEL -1 

The empirical expressions for the scale factor and the scalar field chosen for model-1 are 
given by, 


a = a 0 (£j (11) 

The scale factor (in eqn. 11) produces a constant deceleration parameter, which is given by, 
qo ~ ~ “7 “ “ — 1- Thus a is given by, 


i 

l+Qo 


(13) 


Using equation (11) in (8) we get, 

p = p 0 a 0 - 3(1+r \t/t o r 3 “^ (14) 


Using equations (11), (12) and (14) in equation (6) and writing all parameter values for the 
present epoch (t = t 0 ), one obtains, 



(15) 


Equation (15) is quadratic in n. Two solutions for n are given by, 


5 11 

n+ — - 1 -+ - 

2 2a — 2 



+ 4 p 0 t 0 2 
a 2 (p 0 


(1 - To) 


( 16 ) 


For the present matter dominated era one may take y 0 — 0. 
One thereby obtains its two values as n+ — —1.987, —2.563. 

Using equations (11) and (12) in equation (9), one gets, 
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CO = 
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(17) 
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The time derivative of the BD parameter from equation (17) is given by, 

. _ pot 0 a 0 -3(1+y) {2a(3+3y+n)-4} / t \ 1 -3a(l+y)-na 

(p 0 n 2 a 2 \t 0 J 

Using equations (11) and (12) in equation (10), one gets, 

3 (0.5—1.5y)pot 0 2 a 0 _3( ' 1+ ' / ' ) / t \ 2— 3a(l+y)—na 05 

(a) —-1- ; --- ( — ) - (i)t 

2 <; p 0 na{3a+na-l) \t 0 / 3a+na-l 

Substituting for co in equation (19) from equation (18) one obtains, 

_ 3 1 p 0 t 0 2 a 0 3 ^ 1+y \4-6a-6ay-3nay-na) / t \ 2 ~3a(l+y)-na 

2 2 (p 0 n 2 a 2 (3a+na~l ) \t 0 / 


(18) 


(19) 
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Combining equation (11) with (20), one gets the following expression (eqn. 21) for wasa 
function of the scale factor (a). 


3 1 p 0 t 0 2 a 0 3< ' 1+ ^(4-6a-6ay-3nay-na) 

CO —- 1 --——--- 

2 2 <p 0 n z a z (3a+na-lj 


(-) 

Va 0 / 


2-3a(l+y)-na 
a 


( 21 ) 


Combining equation (12) with (20), one gets the following expression (eqn. 22) for co as a 
function of the scalar field ( cp ). 


_ 3 (l + ) 2-3o;(l+y)-ng 

3 1 p 0 t 0 a 0 (4-6a-6ay-3nay-na) / <p \ na 

2 2 <p 0 n 2 a 2 (3a+na-l) V<p 0 / 


( 22 ) 


MODEL-2 

The empirical expressions for the scale factor and the scalar field chosen for model-2 are 


given by, 


a = a 0 Exp [/? - l)] 

(23) 


(24) 


The scale factor (in eqn. 23) produces a constant deceleration parameter, which is given by, 

9o — ~~7 — ~ 1- Here, the Hubble parameter is // = - = —. Since it is independent of time, 

Ce Cl t q 

it must be equal to the present value of Hubble parameter (. H 0 ). Thus we have, 

/? = H 0 t 0 (25) 

Using equation (23) in (8) we get, 
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P — Po a o 3(1+Y) Exp -3/?(l +y) (“ 1 ) 


(26) 


Using equations (23), (24) and (26) in equation (6) and writing all parameter values for the 
present epoch (t = t 0 ), one obtains, 


© m2 + (?) m + {$“S (1_, '» ) ) = 0 

Equation (27) is quadratic in m. Its two solutions are given by, 


(27) 


m+ = — - ± - l + ^(l- 7 o) 

1 2 - 2d P 2 cp 0 v 


(28) 


For the present matter dominated era one may take y 0 = 0. 
One thereby obtains its two values as m ± — —1.966, —3.034. 

Using equations (23) and (24) in equation (9), one gets, 


O) = 


m 2 /3 2 


3/? 2 (l + m) — 


. 2 —3d+y) 

Potp a o _ 

Vo 


Exp |(— 3/? — 3/?y — m/3) (J- — l)j 


(29) 


The time derivative of the BD parameter from equation (29) is given by, 

* = {(-V - m - mfi) (i - 1)} 

Using equations (23) and (24) in equation (10), one gets, 

3 (l-3y)p 0 to 2 ao _3(1+K) Exp((-3^-3/?y-m/?)(T_ 1 j| to ^ m 

^ 2 2/? z cp 0 m(3+m) 2/? z m(3+m) ^ 

Substituting for d> in equation (31), from equation (30), one obtains, 


3 1 (-6-6y-3my-m)p 0 t 0 2 ao 3( ' 1+y ' ) 

2 2 (p 0 m 2 p 2 (3+m) 


Exp j(— 3 — 3y — m)(3 



(30) 


(31) 


(32) 


Combining equation (23) with (32), one gets the following expression (eqn. 33) for as a 
function of the scale factor (a). 


3 1 (-6-6y-3my-m)p 0 t 0 2 a 0 3 ^ 1+r - > / a \ 3 3y m 

2 2 <p 0 m z /? z (3+m) Va 0 / 


(33) 


Combining equation (24) with (32), one gets the following expression (eqn. 34) for oj as a 
function of the scalar field (< p ). 


3 1 (-6-6y-3my-m)p 0 to 2a o 3( ' 1+,/ ' ) 

2 2 <p 0 m z /? z (3+m) 


-3/?-3/?y-m/J 



(34) 
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3. 


Results and Discussions 


The values of the BD parameter at the present epoch (t = t 0 ), from model-1, for two values 
of n (i.e. n+, as in eqn. 16) are given by, 


3 1 Po t 0 2 a 0 - 3(1 +ro) ' 

2 2 (p 0 n± 2 a 2 (3a+n±a-l ) 


T4-6a-6ayo-3n-j.ayo-n.j-a) 


(35) 


For the present matter dominated era we have taken y 0 — 0. 

Using the values of the parameter n, i.e., n + — —1.987 and n_ = —2.563, 
one gets coo-j. = —1.503, —1.430. 


The values of the BD parameter at the present epoch (t = t 0 ), from model-2, for two values 
of m (i.e. m+, as in equation 28) are given by, 


f \ ? — 3(i+yo) 

_ _ 3 1 (-6-6y-3m ± yo-m ± )poto a 0 _ 

2 2 (p 0 m+ 2 fi 2 (3+m+) 


(36) 


For the present matter dominated era we have taken y 0 — 0. 

Using the values of the parameter m, i.e., m + = —1.966 and m_ = —3.034, 
one gets = —1.518, —1.333. 

According to a study by Banerjee and Pavon, the value of co must be greater than —3/2 [13]. 
This requirement is satisfied if we take n = n_ in model-1 and m = m_ in model-2. 
Therefore, to determine the time dependence of u>, one should use these values of n and m in 
the equations (20) and (32), for model-1 and model-2 respectively. 

Using all parameters values and taking n — n_, the expression of the BD parameter m (eqn. 
20), based on model-1, takes the following form. 

a) = -1.5 + (0.071 - 0.073y) (^J (37) 


From equation (37), the time derivative of a) is given by, 


du> 

dt 


(0.073 — 0.548y+0.487y 2 )(A) 


1.029-6.667# 


t 


(38) 


From equation (38), the value of — at the present epoch (t = t 0 ) is obtained as, 


(dcj\ _ 0.073-0.548y„ + 0.487y o 2 
\dt) t=to t 0 


(39) 


Using all parameters values and taking m = m_, the expression for the BD parameter co 
(eqn. 32), based on model-2, takes the following form. 


a) = -1.5 + (0.168 - 0.175y)£xp [(0.035 


S.MSy)^- 1 )] 


(40) 


From equation (40), the time derivative of a) is given by, 
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(41) 


da) (o. 006-0.518Y+0.533y 2 ) Exp[(0. 035-3. 045y)(^—l)] 

dt to 


From equation (41), the value of — at the present epoch (t — t Q ) is obtained as, 


_ 0.006-0.518y 0 +0.533y 0 2 
\dtJ t=to t 0 


(42) 


The expressions of (39) and (42), based on model-1 and model-2 respectively, are identical in 
form. Here, t 0 denotes the age of the universe, with a very large value of 1.4 x 10 10 years or 


4.42 x 10 17 seconds. Due to this extremely large value, 



would be very small. The 


sign of 



is determined by the value of y of the cosmic fluid at the present epoch. 


This value has been taken to be zero for the present era of matter dominated universe 
(pressureless dust) by Sahoo and many researchers [17]. Setting the right hand sides of 
equations (39) and (42) equal to each other, one gets, 


0.046y 0 2 + 0.030y 0 - 0.067 = 0 


(43) 


Equation (43) is quadratic in y 0 , leading to the following solutions. 


(7o)± 


— 0.030±V(0.030) 2 —4X0.046X(—0.067) 
2X0.046 


0.924,-1.576 


(44) 


The negative value of y 0 , in equation (44), is consistent with recent observations [19, 20]. 

The results of equation (44) have their origins in equations (38) and (41). We took y 0 — 0 to 
obtain the values of n and m , used for evaluating the constants in the equations (38) and (41). 
We should therefore discuss a more generalized method to estimate the value of the EoS 
(equation of state) parameter y 0 . 

According to this new method, one may determine ( —) for model-1 and model-2, from 

^ dt 't=t 0 

equations (20) and (32) respectively. Since they represent the time dependence of the same 
entity (a>) at the present epoch, their values are likely to be equal to each other. Setting them 
equal to each other and taking m — n, one gets, 

Ay 0 2 +By 0 + C = 0 (45) 

where, 


A = 18a 2 fin 2 — 54a 3 /?n 2 + 54a 2 /? 2 n 2 + 9 a 2 fin 3 — 45 a 3 /?n 3 + 45 a 2 (3 2 n 3 — 9 a 3 fin 4 
H- 9a 2 /? 2 n 4 

B = 36a 2 (In 2 — 108 a 3 /?n 2 — 72 a/3 2 n 2 + 108 a 2 /3 2 n 2 + 18a 2 /3n 3 — 90a 3 /3n 3 — 42 a/3 2 n 3 
H- 90 a 2 /3 2 n 3 + 3a 2 /?n 4 — 2 7a 3 (In 4 — 6 a/3 2 n 4 + 27a 2 /? 2 n 4 — 3a 3 /?n 5 
+ 3 a 2 /3 2 n 5 

C = 18a 2 fin 2 — 54a 3 /?n 2 + 24 /? 2 n 2 — 72 a/3 2 n 2 + 54 a 2 /3 2 n 2 + 9a 2 /?n 3 — 45 a 3 /?n 3 + 8 /? 2 n 3 
— 42 a/? 2 n 3 + 45 a 2 /3 2 n 3 + a 2 /3n 4 — 12a 3 /?n 4 — 6a/? 2 n 4 + 12 a 2 /3 2 n 4 — a 3 (3n s 
+ a 2 /3 2 n s 
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Equation (45) is quadratic in y 0 . 

Here, the values of a and /? can be obtained from equations (13) and (25) respectively. 

The reason for taking n — m is due to the closeness of their values, obtained from equations 
(16) and (28). Since the expressions of cp in equations (12) and (24) are not the exact 
solutions of the field equations, the values of n and m can also be other than those given by 
equations (16) and (28) respectively. The only restriction imposed upon their choice of values 
is the requirement that —1.5 < <n 0 < 0, as obtained by Banerjee and Pavon [13]. We have 
only negative values n and m in equations (16) and (28). Negative values of n and m indicate 
that the scalar field (<p) decreases with time, as evident from their defining equations (eqns. 
12, 24). Consequently, the gravitational constant G = l/cp increases with time. There are 
other recent studies showing a decrease of (p (or an increase of G) with time [17, 18]. 


The values of y 0 obtained from quation (45) are given by, 


Y o 


± 


-B±VB 2 -4i4C 

2A 


(46) 


For — 2 < n < 0, one finds from equation (46) that y 0 _ vary approximately in the range of 
— 1.9 < y 0 _ — 0.6, which is consistent with the recent experimental findings regarding the 
EoS parameter [19, 20]. The value of y 0+ is positive for this range of n values. These positive 
values are not at all consistent with the range of values for y 0 obtained from scientific studies 
based on recent astrophysical observations [19, 20] 


4. Conclusions 

In the present study we have used scale factors that generate time independent deceleration 
parameters. As per scientific studies based on recent observations, the deceleration parameter 
has changed its sign with time, from positive to negative, indicating a change of phase of the 
cosmic acceleration from deceleration to acceleration [18]. Our scale factor for model-1 is 
identical in form to that used by Sahoo et al. [17], although our method of deriving the time 
dependent expression for oo is different. These scale factors can be used for the simplicity of 
mathematical analysis in a theoretical study. They are required to satisfy a condition that the 
constant deceleration parameter obtained from them is negative and q 0 > —1. The time 
dependence of a), for different eras of the expanding universe (characterized by different 
values of the EoS parameter y), is found to be qualitatively the same for these two models. 
The choice of the scalar field in both models is such that it has the same functional 
dependence on the corresponding scale factor. According to generalized Brans-Dicke theory, 
the BD parameter (a>) is a function of the scalar field (cp). In both models, the choices of 
empirical relations regarding the scale factors and the scalar fields have allowed us to 
formulate expressions of oo as functions of cp. Using the results of these two models together, 
the range of values for y 0 (EoS parameter at the present epoch) has been estimated. One 
could also have taken n A m for a more generalized form, leading to a larger mathematical 
expression which might be untenable. An improvement over this method of determination of 
the EoS parameter can be made by using scale factors that are obtained as solutions of the 
field equations for different logical choices of the scalar field. 
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